Avoiding the polarization catastrophe in LaAlC>3 overlayers on SrTiO3(001) 

through a polar distortion 



Rossitza Pentcheva a and Warren E. Pickett fc 
a Department of Earth and Environmental Sciences, 
University of Munich, Theresienstr. 1^1, 80333 Munich, Germany and 
b Department of Physics, University of California, 
Davis, One Shields Avenue, Davis, CA 95616, U.S.A. 
(Dated: December 2, 2008) 

A pronounced uniform polar distortion extending over several unit cells enables thin LaA103 
overlayers on SrTiO3(001) to counteract the charge dipole and thereby neutralize the "polarization 
catastrophe" that is suggested by simple ion-counting. This unanticipated mechanism, obtained 
from density functional theory calculations, allows several unit cells of the LaAlOs overlayer to 
remain insulating (hence, fully ionic). The band gap of the system, defined by occupied O 2p states 
at the surface and unoccupied Ti 3d states at the interface in some cases ~20 A distant, decreases 
with increasing thickness of the LaA103-film before an insulator-to-metal transition and a crossover 
to an electronic reconstruction occurs at around five monolayers of LaAlC>3. 

PACS numbers: 77.22.Ej,71.28.+d,73.20.Hb,75.70.Cn 



Heterostructures containing polar (e.g. LaA103 
(LAO)) and nonpolar (e.g. SrTi0 3 (STO)) oxides 
raise issues similar to semiconductor heterostruc- 
tures, e.g. Ge/GaAs If the ionic charges are 
maintained, then adding LAO layers (with alternat- 
ing positively charged LaO and negatively charged 
A102-layers) onto an STO substrate is expected to 
lead to an ever-increasing dipole. The accompa- 
nying potential shift across the LAO slab should 
eventually cause a "polarization catastrophe" @, Q . 
To avoid this, surfaces are supposed to reconstruct 
(atomically, e.g. via defects) or facet Q. However, 
in oxides containing transition metal ions further 
compensation mechanisms are available. These can 
lead to electronic phenomena that are unanticipated 
from the properties of the bulk constituents and can 
thus invoke new functionality. Although LAO and 
STO are two conventional nonmagnetic band insu- 
lators, a high conductivity H as well as indications 
for magnetism Q and superconductivity Q were 
measured at the n-type LaO/Ti02 interface (IF). 
The initially reported high carrier density was sub- 
sequently found to be sensitive to film growth con- 
ditions [2, H, H, [H, [l(| EH ■ Suppressing oxygen va- 
cancies leads to an increase of the sheet resistance 
by several orders of magnitude [f| [l2j . 

Recently, a thickness dependent switching be- 
tween insulating and conducting behavior [13[ was 
reported around four monolayers (MLs) of LAO on 
STO(OOl). Thiol et al. [3 found that the insulator- 
to-metal transition can also be induced dynamically 
by an external electric field. Writing and erasing of 
nanoscale conducting regions in LAO overlayers on 
SrTiO3(001) was demonstrated [1J| with an applied 
local electric field (by atomic force microscope). In 
contrast to the abrupt transition in LAO films on 



STO(OOl), Huijben and collaborators [15| observed 
a smoother increase in carrier density in coupled 
p- (hole-doped) and n-type (electron-doped) inter- 
faces where the LAO film was covered by SrTi03. 
These results lend additional urgency to understand 
the insulator-to-metal switching in oxide nanostruc- 
tures. 

The seeming violation of charge neutrality has 
become a central consideration in trying to under- 
stand the unexpected behavior that emerges at the 
interface. So far, theoretical efforts have concen- 
trated mainly on LAO/STO superlattices [H, [13, 
H [H |M E], |H [H. For these, density functional 
theory calculations (DFT) suggest that the charge 
mismatch at the electron doped n-type interface is 
accommodated by partial occupation of the Ti 3d 
band. At low temperatures a disproportionation on 
the Ti sublattice in the interface layer is predicted 
within LDA(GGA)+U with Ti 3+ and Ti 4+ ar- 
ranged in a checkerboard manner 

At first glance the transition from insulating to 
conducting behavior found between three and four 
MLs on SrTiO 3 (001) does not fit into the picture 
obtained from the GGA+U calculations for infinitely 
extended superlattices [lfl H^|. In the latter the 
charge mismatch is accommodated immediately at 
the IF and largely independent of the thickness of 
LAO and STO. 

Here we shed light on the nature of the electronic 
state at the interface and the origin of thickness 
dependent metal-to-insulator transition. Based on 
DFT-calculations, we demonstrate that an interface 
between a thin LAO film on a STO (001) substrate 
is fundamentally different from an IF in a periodic 
superlattice. In a thin overlayer there are two polar 
discontinuities, one at the IF and another at the sur- 



face. As we will show in this paper, the proximity 
to the surface and the associated atomic relaxation 
has a profound effect on the properties and compen- 
sation mechanism at the interface. 

To address these issues, we study here the struc- 
tural and electronic properties of 1-5 ML of LAO on 
a STO(001)-substrate, and contrast them with the 
corresponding behavior at LAO/STO superlattices. 
The DFT calculations are performed with the all- 
electron full-potential augmented plane waves (FP 



APW) method in the WIEN2k-implementation 25 1 
and generalized gradient approximation (GGA) 26 1 
of the exchange-correlation potential. We have used 
the lateral lattice constant of the SrTiC>3-substrate 
obtained from GGA (3.92A) because it is in closer 
agreement with the experimental value (3.905A). 
Thus the epitaxial LaA103-film (3.79A) is subject to 
tensile strain on SrTiOs due to a lattice mismatch of 
3%. The GGA-gap for bulk SrTi0 3 is 2.0 eV (exper- 
imental value 3.2 eV) separating filled O 2p bands 
and unfilled Ti 3d bands. For LaA103 we obtain a 
band gap within GGA of 3.7 eV (experimental value 
5.6 eV) between filled O 2p bands hybridized with 
Al p-bands and unfilled conduction bands comprised 
of La 5d and Al 3s, 3p states. 

The thin (1-5 ML) LAO films on a STO(OOl) sub- 
strate were modeled in the supercell geometry with 
two inversion symmetric surfaces to avoid spurious 
electric fields. To minimize the interaction between 
neighboring surfaces, the slabs are separated in z- 
direction by 10-12A of vacuum. Because XRD [27j 
and LEED [28| measurements give no indication for 
a superstructure, we have considered only defect-free 
unreconstructed surfaces. However, we have allowed 
for more complex symmetry breaking and electronic 
reconstruction by using a c(2 x 2)-lateral periodicity. 
For this cell 15 fc-points in the irreducible part of the 
Brillouin zone were used. 

Polar distortions. A full structural optimization 
was performed within the tetragonal unit cell using 
GGA. The relaxed structures are pictured in Fig. Q] 
and the atomic displacements which are exclusively 
along the c-axis are plotted in Fig. [2] 

An unexpected pattern is uncovered: (i) The sur- 
face AIO2 layer shows an inward relaxation of both 
Al (0.12-0. 19A) and O (0.11-14A) atoms by simi- 
lar amounts, resulting in only a small buckling of 
0.01-0.04A. (ii) The most striking result is a strong 
buckling in all LaO layers with a uniform polar-type 
distortion of 0.20-0.36 A. This relaxation is domi- 
nated by the outward movement of La by up to 0.26 
A. (Hi) The subsurface AIO2 layers are also buck- 
led by approx. 0.15 A: Al relaxes outward by 0.01- 
0.08 A, O inward by - 0.10 A. (iv) The Ti and 
O-ions at the interface show only small outward re- 
laxation (less than 0.06 A). This situation of large 
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FIG. 1: Side view of the relaxed structures of 1, 3 and 4 
ML LAO on STO(001) showing the ferroelectric distor- 
tion. The oxygen ions are marked by light grey spheres, 
while the St-, Ti-, La- and Al-cations are shown by green, 
blue, purple and orange spheres. 



polar distortions throughout the LAO slab with neg- 
ligible distortion in STO contrasts strongly the be- 
havior obtained for LaA10 3 (or LaTi0 3 )/SrTi0 3 
superlattices where ferroelectric-like plane-buckling 
is predominantly in the Ti02 layer at the inter- 
face 0, E El M, US HI and amounted to 0.17A 
for a 5LAO/5STO heterostructure 

Although not explicitly discussed, previous VASP 
calculations showed a similar relaxation pattern [l4| . 
In contrast, another WIEN2k study reported only 
small relaxations [3l[ also finding metallic behavior 
for all LAO thicknesses. Vonk et al.^j^ fitted struc- 
tural models to surface x-ray diffraction (SXRD) 
data for a LAO overlayer on STO(001) with a total 
coverage of 0.5 MLs. The best fit infers a contracted 
TiOg octahedron at the interface and an outward re- 
laxation of La in the LaO layer. A direct comparison 
to our results is not possible because of the different 
coverage and since the incomplete LAO layer allows 
for lateral relaxations of strain. 

Compensating electronic and ionic dipoles. The 
implied dipole moment of a 4ML LaA103-film is 
Vbare ~ -0.5e x (3.9A) x 4cells = -7.8 e°A which 
translates into a bare potential shift across the four 
LAO layer slab of 80-90 eV. This electric field is 
screened by electronic plus ionic polarization. Using 
the bulk LAO dielectric constant e—2A, this becomes 



approximately AV = Aire 



3.5 eV, which 



e(3.9A) 2 

can be sustained by the 5.6 eV gap of LAO. 

The dipole shift due to the polar distortion can be 
estimated from the displacements, using the formal 
ionic charges (33[. The layer resolved dipoles, pic- 
tured in Fig. [2l show a strikingly uniform behavior 
in the inner LAO layers (a strong dipole moment of 
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FIG. 2: Vertical displacements of ions Az with respect 
to the bulk positions in A for 1-5 MLs of LAO on 
STO(OOl) obtained within GGA. a) and b) displays the 
relaxations in the AO and BO2 layers, respectively with 
cation/oxygen relaxations marked by full black/open red 
symbols. The corresponding layer- resolved dipole mo- 
ments are displayed in (c). The x axis shows the distance 
from the interface (I) Ti02-layer. 

~ 0.55 and ~ 0.75 eA in the AIO2 and LaO layers, 
respectively) and an almost vanishing dipole in the 
interface TiC>2 and surface AIO2 layers. The total 
ionic dipole e.g. for the 4ML LAO film is: 

V = ZiAzi » 4.8 eA (1) 

This rough estimate gives a value that is more than 
60% of Vtare and of opposite sign. The large polar 
distortion has a significant impact on the electronic 
behavior of the system (LAO film + STO substrate). 

Electronic properties of overlayer + interface. 
The total density of states as a function of the LAO 
thickness for 1-5 ML thick LAO films on STO(OOl) 
is shown in Fig. [3] For the ideal atomic positions 
all systems are metallic due to an overlap of the sur- 
face O 2p bands with the Ti 3d bands at the in- 
terface (cf. Fig. [4]). The strong lattice polariza- 
tion leads to insulating behavior for N = 1 — 4 ML. 
With increasing thickness of the LAO-film the gap 
decreases by approximately 0.4 eV per added LAO- 
layer starting from 1.7 eV for 1ML LAO and finally 
closes for 5 ML LAO/STO(001). While the DFT 
calculations correctly predict an insulator-to-metal 
transition, the critical thickness is affected by the 
band gap underestimation within GGA and is ex- 
pected to occur beyond an LAO thickness of 6 MLs 



FIG. 3: Density of states for the ideal (dashed line) and 
relaxed (solid line, grey filling) structure of 1-5 ML LAO 
on STO(001). Relaxation opens a band gap, but its size 
decreses with each added LAO layer. 



for a defect-free surface. Defects and/or adsorbates 
which are not considered in the theoretical model 
may be responsible for the lower critical thickness 
found in experiment [l3T |. 

In order to gain more insight into the type of elec- 
tronic states around the Fermi level and their spatial 
distribution, we have plotted in Fig. 0] the layer- 
resolved density of states for 5ML LAO/STO(001) 
with ideal and relaxed atomic positions. A nearly 
rigid upward shift of the O 2p-bands towards the 
Fermi level is observed as they approach the sur- 
face. For the relaxed configuration the top of the 
valence band is determined by the O 2p band in the 
surface AIO2 layer, while the bottom of the conduc- 
tion band is fixed by the unoccupied Ti 3d band at 
the interface. Thus, the closing of the gap is due to 
the raising of the O 2p valence band states in the 
surface LAO layer as the LAO slab gets thicker. 

For an isolated IF (in LAO/STO superlattices), 
the "extra" 0.5e~ per IF Ti ion must be accommo- 
dated, it cannot go into the LAO layers due to the 
large band gap, and it remains at the IF and charge- 
orders (within GGA+U) rather than occupy itiner- 
ant Ti 3d states within the STO layers [l&Slj]. In 
contrast, in thin LAO overlayers, it does not go into 
the Ti 3d states anywhere; the lattice screening pre- 
serves the formal electronic charges corresponding 
to closed shells, therefore including correlation ef- 
fects within LDA+U would not influence the result. 
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FIG. 4: Layer resolved density of states of 5 ML LAO 
on STO(OOl) with ideal (dashed line ) and relaxed (grey 
shaded area) coordinates. The DOS for the ideal posi- 
tions was shifted by 0.5 eV to align with the conduction 
band of the system with the relaxed atomic positions. 
The Fermi level of the unrelaxed system is marked with 
a dashed line. 

Up to a critical LAO thickness of 5ML the 'extra' 0.5 
electron at the interface appears to have been shifted 



to the surface and there is no need for an electronic 
(or atomic) reconstruction as observed for the su- 
perlattice. However, beyond a critical coverage, the 
system is expected to become metallic. 

In summary, strong lattice polarization compen- 
sates the dipolar electric field in thin LAO-overlayers 
and sustains the insulating behavior up to 5 MLs 
LAO. The effect of this polar distortion is crucial: 
Without it, the system is metallic for any coverage 
of LaA103. The band gap in the relaxed systems 
formed between the 02p states in the surface layer 
and the Ti 3d states at the interface decreases grad- 
ually and collapses at 5 ML LAO/STO(001) giving 
rise to an electronic reconstruction as the system ap- 
proaches the limit of the isolated interface. 

The structural distortion we observe here dif- 
fers from issues discussed for ultrathin ferroelec- 
tric films, e.g. the existence of a critical thickness 
below which ferroelectric displacements are sup- 
pressed [lil |H, [H, H3]- In contrast, LaA103 is a 
wide band gap insulator with a low dielectric con- 
stant and no ferroelectric instability. In this case 
the lattice distortion is driven by the polar nature 
of the surface. Further examples that nanoscale ob- 
jects can sustain polarity are recently emerging (see 
[381 ] and references therein). In the context of tai- 
loring materials properties for device applications, 
finite size effects in ultrathin oxide films prove to be 
an exciting area for future research. 
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